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An orthogonal surface-oriented 
coordinate system has been developed 
for three-dimensional flows where the 
computational domain normal to the 
surface is small. With this restric- 
tion the coordinate system requires 
orthogonality only at the body surface. 
The coordinate system is as follows; 
one coordinate measures distance 
normal to the surface while the other 
two coordinates are defined by an 
orthogonal mesh on the surface. One 
coordinate is formed by the inter- 
section of the body surface and the 
meridional planes as illustrated in 
Figure 1 and gives the 9 = constant 
lines. The other coordinate which 
is nondimensionalized with a character- 
istic length of the body geometry, 
measures the distance along the body 
surface when 0=0. This coordinate 
system has been utilized in boundary 
layer flows''^ and for the hypersonic 
viscous shock- layer problem.^ 

Two methods have been developed 
for generating the surface coordinates. 
The first method uses the orthogonal 
condition in finite-difference form to 
determine the surface coordinates with 
the metric coefficients and curvature 
of the coordinate lines calculated 
numerically. The second method obtains 
analytical expressions for the metric 
coefficients and for the curvature of 
the coordinate lines. Both methods 
assume the body surface is defined in 
terms of a cylindrical coordinate 
system where r = r(x,fl). The surface 
inclinations and as illustrated 
in Figure 2 are determined from 

tan = (If and r tan = - (|f)^ 
and are known quantities. 


PLANE OF I 



Figure 1. Surface Coordinate System. 



Figure 2. Angles <!>, and Defined in 
the Cylindrical Coordinate 
System. 
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In the numerical method,*'^ the 
orthogonal condition for the surface 
coordinates results in the relation 


1 



2 

- 5 —^ = 2TTr cos <fl- 

d'i) 1 


tan 



dx = A de (along C = constant) 

where 

2 

A = r tan <ti^ tan 412/(1 + tan 4j) 

The equation of the surface provides 
the relation 

dr = tan 4^ dx - r tan 42 <36 

The surface coordinate C is determined 
numerically from the foregoing 
equations by assuming a value of d6 
and marching away from 6 = 0 to deter- 
mine the values of x and r. In 
addition the metric coefficients are 
determined numerically from 

h^ = ds/dS 
h = dt/do 

CjJ 

where 


This equation is integrated along 
C = constant lines on the surface 
from the initial condition hr = 

1 at 0 = 0 . The substitution of 
foregoing equations into the 
equations for and give 
analytical expressions for the 
curvature of the coordinate lines. 
In evaluating these relations, 
the variations of x and 6 along 
the ^ = constant coordinate must 
be known. 

A sphere at angle of attack 
as shown in Figure 3 is used 
to illustrate the computation of 
the surface coordinates with both 
methods. The surface coordinates 
on the sphere as viewed from the 
side are illustrated in Figure 4. 
The C = constant lines result from 
planes intersecting the sphere 
with these planes passing through 
the line which is normal to the 
plane of symmetry and is located 
at 


U) = 

6/2TT 


ds^ = 

dx^ + 

dr2 

dt^ = 

ds^ + 

r^dO 


The curvature of the coordinate lines 
is determined from 

„ _ 1 for 0 ) = constant 

C “ h^h^ 

„ 1 "to for C = constant 

u “ h.h 3C 
C to 


with the derivatives replaced with mid- 
point difference relations. 

In the second method^, an analyt- 
ical expression is developed for h 
as follows: “ 

h = 2trr (1 + cos^4i tan^4^)^^^ 
to I i 

A differential equation results for 
the other metric coefficient as follows: 


x/a = V 1 - (b/a) ^ 

y/a = (x/a) ^/(b/a) 

The metric coefficients for this 
coordinate system are given in 
Figures 5 and 6 with good agree- 
ment between the two methods. The 
curvature of the coordinate lines 
is given in Figures 7 and 8 . It 
is noteworthy that K, is independent 
of C. The differences evidenced 
in Figure 8 can be partially 
attributed to the numerical evalua- 
tion of being at one-half mesh 
space locations away from the C 
indicated. 

The numerical method of gene- 
rating the orthogonal surface 
coordinates has been applied to 
ellipsoids, paraboloids and elliptic- 
paraboloids. The coordinates on an 
ellipsoid are illustrated in Figure 9. 
The second method or analytical 
approach has only been developed for 
the sphere. 







Figure 7. Geodesic Curvature of 
Lines of Constant 6. 
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Figure 9. Surface Coordinates on 
Ellipsoid (b/a = 1/4) 


Figure 8. Geodesic Curvature of 
Lines of Constant C- 
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